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Flows of suspensions are often affected by wall slip, that is the fluid velocity Vf in the vicinity of a 
boundary differs from the wall velocity Vw due to the presence of a lubrication layer. While the slip 
velocity Vs = \vf — v.ui\ robustly scales linearly with the stress cr at the wall in dilute suspensions, 
there is no consensus regarding denser suspensions that are sheared in the bulk, for which slip 
velocities have been reported to scale as a Us oc with exponents p inconsistently ranging between 
0 and 2. Here we focus on a suspension of soft thermoresponsive particles and show that Vs actually 
scales as a power law of the viscous stress a — a^, where CTc denotes the yield stress of the bulk 
material. By tuning the temperature across the jamming transition, we further demonstrate that 
this scaling holds true over a large range of packing fractions (j> on both sides of the jamming point 
and that the exponent p increases continuously with cj), from p = 1 in the case of dilute suspensions 
to p = 2 for jammed assemblies. These results allow us to successfully revisit inconsistent data from 
the literature and paves the way for a continuous description of wall slip above and below jamming. 

PACS numbers: 83.50.Rp, 83.80.Hj, 83.60.La, 47.57.E- 


I. INTRODUCTION 

Wall slip is ubiquitous in flows of complex fluids and 
plays a key role in numerous situations ranging from tur¬ 
bulent drag reduction to extrusion processes [T]. It de¬ 
notes a situation where the fluid velocity Vf in the vicinity 
of a boundary differs from the velocity Vm of the bound¬ 
ary itself, and often originates from a chemical or physical 
mismatch between the constituents of the fluid and that 
of the wall. Such apparent discontinuity in velocity at 
a boundary reflects the existence of a thin and highly 
sheared region adjacent to the wall where the viscosity 
is much smaller than in the bulk material [5]. Although 
historically described as a mere artifact, wall slip appears 
as an intrinsic feature of complex fluids that has been re¬ 
ported in polymeric fluids ma as well as in soft glassy 
materials such as colloidal suspensions liH, foams i- 
[TO] . emulsions [IIHI3], microgels [mill], attractive gels 
M: etc. Moreover, wall slip is often non-trivially cou¬ 
pled to the bulk dynamics of complex fluids, especially 
during the yielding transition where it may even govern 
the nature of the subsequent steady-state flow [nHTi]. 

Physical insights into wall slip are generally gained 
from the scaling of the slip velocity Vg = \vf — Vw\ with 
the stress at the wall a. More specifically for particulate 
suspensions below jamming, the particle deformability 
does not play a key role and slip velocities scale roughly 
linearly with a [iiiiiDHia. In jammed particulate as¬ 
semblies, however, the situation is much more subtle as 
the material exhibits a yield stress CTc below which it 
cannot flow in the bulk. The main goal of the present 
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Rapid Communication is to clarify the slip behavior of 
a jammed suspension of soft particles driven above the 
yield stress. Indeed in this case Vg has been reported 
to increase as a quadratic function of a in concentrated 
emulsions and microgels [mils] while other experiments 
on similar samples rather led to a linear scaling of Vg with 

a [HIMIIIS]. 

Therefore despite a large body of work on wall slip, the 
key parameters controlling the scaling of Vg in a dense, 
yielded suspension are still unknown. One of these pa¬ 
rameters is the packing fraction which effect on slip has 
been scarcely investigated. Here we focus on a suspension 
of thermoresponsive poly(Af-isopropylacrylamide) parti¬ 
cles that allows us to explore the slip behavior of a sus¬ 
pension at various packing fractions in a single set of ex¬ 
periments, from a dilute system up to a jammed assem¬ 
bly of soft particles. Extensive rheological experiments 
coupled to velocimetry demonstrate that the slip veloc¬ 
ity scales as a power-law on both sides of the jamming 
transition. For low packing fractions we recover Vg oc 
with p ~ I, while above jamming and for stresses larger 
than the yield stress, Vg follows a power law of the viscous 
stress, i.e. Vg ot {a — CTc)^. Furthermore the exponent p 
increases continuously from p = 1 to p = 2 as the packing 
fraction is increased. This is confirmed by revisiting for¬ 
merly published data on emulsions and microgels. Our 
results thus explain long-standing discrepancies in the lit¬ 
erature and bridge the gap between the slip behavior of 
dilute suspensions and soft jammed materials. 


II. EXPERIMENTAL 

pNIPAM microgels.- Samples consist in poly(A- 
isopropylacrylamide) or pNIPAm microgel particles 
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FIG. 1: (color online) (a) Average diameter of the pNIPAm 
particles vs temperature T as measured by dynamic light scat¬ 
tering. (b) Elastic (G', ■) and viscous {G”, □) moduli mea¬ 
sured by small amplitude oscillatory shear (7 = 1 % and 
uj = 1 Hz). The data correspond to the average of an increas¬ 
ing and decreasing ramp of temperature of duration 25 min 
each, (c) and (d) Flow curves, shear stress a versus shear 
rate 7, recorded at different temperatures: T = 17, 20, 26 
and 34° C from top to bottom. Each flow curve is obtained by 
sweeping down 7 from 500 to 5.10“® s“^ with a waiting time 
of 10 s per point. The red curves are the best fits of the data 
by the Herschel-Bulkley model with = 45.7T0.2, 31.3±0.4 
and 8.9 ±0.1 Pa from top to bottom in (c) and by a power 
law fit, a = with n = 0.74 and 77 = 0.0245 Pa.s®'^^ in 
(d). 


synthesized by aqueous free-radical polymerisa¬ 
tion from N, A^'-methylenebisacrylamide and N- 
isopropylacrylamide (Sigma Aldrich) following a recipe 
previously described in and more recently 

applied in [29] . The reaction is performed in a three- 
neck round-bottom flask (500 mL), equipped with a 
magnetic stir bar, a reflux condenser, a thermometer 
and an argon inlet. 7V-isopropylacrylamide (NIPAm) 
monomers at an initial concentration of 70 mM and 
N, A^'-methylenebisacrylamide (BIS) at 2.5 mol% are 
dissolved in 98 mL of water. The solution is purified 
through a 0.2 /rm membrane filter to remove residual 
particulate matter, before being heated up to 70°C with 
argon thoroughly bubbling. Free radical polymerization 
is initiated after 1 h, with potassium persulfate (KPS, 
2.5 mM) dissolved in 2 mL of water. The solution 
initially transparent becomes progressively turbid due 
to the polymerization process which is carried out 
during 6 h in the presence of argon and under stir¬ 
ring. To eliminate chemical residues and water-soluble 
oligomers, the microgels are then purified by at least 
5 centrifugation-redispersion cycles (under 21.10® g for 
1 h, where g denotes the acceleration of gravity). For 


each cycle, the supernatant is removed and its surface 
tension is measured by the pendant drop method. The 
purification process is repeated until the surface tension 
of the supernatant is comparable to that of pure water, 
i.e. above 70 mN/m. After the last centrifugation 
step, the supernatant is totally removed. The microgel 
concentration is quantified by measuring the dried mass 
of polymer in the paste and estimated to be equal to 
5.2% wt. The polydispersity index of the sample was 
estimated through Dynamic Light Scattering to be con¬ 
stant and about 0.06 in the jammed state (T < 32°C). 
Such a polydispersity should be enough to ensure that 
no crystallization is taking place. Finally, the microgel 
suspensions are seeded with acoustic tracers incorpo¬ 
rated at 35°C, namely 1% wt polystyrene microbeads 
(Dynoseeds TS, 20 diameter, density 1.05), that 
were checked to have a negligible influence on the system 
rheology. 

Velocimetry setup.- Simple shear experiments are per¬ 
formed in a polished Plexiglas Couette cell (typical 
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FIG. 2: (color online) Steady-state velocity profiles v{r) where 
r is the distance to the moving wall: (a)-(c) at T = 17° G for 
different shear rates 7 = 0 . 1 , 10 and 100 s“® from top to 
bottom and (d)-(f) at 7 = 0.1 s“^ for different temperatures 
T = 20, 26 and 34°G from top to bottom. Blue (resp. red) 
lines are the best linear fits of the velocity profiles over 50 to 
150 jj,m from the rotor (resp. stator). Extrapolating each fit 
to the wall gives the fluid velocity u/ at both boundaries. In 
all cases the upper limit of the vertical axis corresponds to 
the rotor velocity. 
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FIG. 3: (color online) (a)-(d) Steady-state slip velocities Vs measured at the rotor (□) and at the stator (■) for different 
temperatures T = 17, 20, 26 and 34°C. Gray lines are the best power-law fits of the data, Vs = B{a — o^Y, where is fixed to 
the yield stress inferred from the Herschel-Bulkley fit of the flow curve (Fig.j^ and B and p are free parameters, (e) Exponent 
p vs temperature T obtained from the previous fitting procedure (•) and by using the value of Cc that minimizes the of 
linear fits of InUs vs ln((j — Uc) when Uc is varied (o). In the former case error bars are computed from the uncertainty on 
(7c, whereas in the latter case error bars correspond to the standard deviation associated with the least-square method. The 
horizontal dashed lines emphasize p — 1 and 2. 


roughness 15 nm, height H = 28 mm, rotating inner 
cylinder of radius R = 23.96 mm, fixed outer cylinder 
of radius 25 mm, gap e = 1.04 mm) equipped with a 
homemade lid to minimize evaporation. Temperature is 
controlled within 0 . 1 °C by a water circulation around 
the Couette cell. Rheological data are recorded with 
a stress-controlled rheometer (MGR 301, Anton Paar) 
while the azimuthal velocity v is measured simultane¬ 
ously as a function of the radial distance r to the rotor 
at about 15 mm from the cell bottom and with a spa¬ 
tial resolution of 40 /im by means of ultrasonic velocime- 
try [30]. For each temperature, velocity measurements 
require the knowledge of (i) the sound speed Cg in the 
sample, which is deduced from ultrasonic time-of-flight 
measurements using a standard transmission setup (Ta¬ 
ble |T]) and (ii) the position of the stator-fluid interface 
and the incidence angle of the ultrasonic beam, that we 
both determine by calibration in an aqueous dispersion 
of 1% wt polystyrene microbeads [30] . 

III. RESULTS 

Properties of the thermoresponsive eolloids.- The indi¬ 
vidual size of pNIPAm microgel particles is known to be 
continuously tunable with temperature, which can thus 
be used to vary the packing fraction of the suspension 
[SIII32]- Above 34° C the microgels are collapsed and the 


sample behaves as a dilute, liquid-like suspension of hard 
particles with a shear-thinning behaviour [Fig. [ 3 a,d)]. 
As the temperature is decreased the particles swell and 
their diameter d increases from 300 to 600 nm over a 
temperature span of about 4°C around the lower criti¬ 
cal solution temperature (LOST) Tic — 32°C [Fig. [^a)]. 
As a result the sample jams into a solidlike material for 
T < Tic [Fig. [^b)] [33l[34]- The flow curves, shear stress 
a versus shear rate 7 , show that the jammed suspensions 
behave as yield stress fluids under shear and are well de¬ 
scribed by the Herschel-Bulkley model, cr = CTc + A 7 ", 
with a yield stress CTc that steeply decreases as T ap¬ 
proaches Tic [Fig. [^c) and Table |l]. 

Slip velocity measurements. - Once the sample is loaded 
into the shear cell, we wait for its viscoelastic moduli to 
reach a plateau. We then impose various shear rates 7 
and record the velocity profiles together with the stress 
response before moving to another temperature and re¬ 
peating the same protocol. At a given temperature, two 
consecutive shear start-up experiments are separated by 
a preshear at 50 s“^ during 40 s to erase the sample his¬ 
tory, followed by a 3 min rest period during which we 
monitor the evolution of G' and G" through small am¬ 
plitude oscillatory shear. Quantitative agreement in the 
viscoelastic moduli measured between successive shear 
experiments indicates that the sample neither ages nor 
suffers from evaporation issues. Steady-state velocity 
profiles obtained at T = 17° C for different shear rates 
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T (°C) 

Cs (m.s ^) 

ac (Pa) 

n 

A (Pa.s") 

(Pa) 

17 

1533±4 

45.7±0.2 

0.44±0.02 

15.5 ±0.2 

45.8±0.1 

20 

1536±4 

31.3±0.4 

0.49±0.02 

8.4 ±0.2 

30.8±0.1 

26 

1541±5 

8.9±0.1 

0.54±0.02 

2.33±0.07 

8.7±0.05 

34 

1548±5 

0 

0.74±0.02 

0.23±0.01 

0 


TABLE I: Sound speed and rheological parameters of the pNI- 
PAm microgel suspensions at different temperatures T. The 
last column shows the value of Cc inferred from the least- 
square procedure detailed in the main text. 

are reported in Fig.[^a-c). The choice of smooth bound¬ 
ary conditions favors wall slip over a large range of shear 
stresses as evidenced by the difference between the ex¬ 
trapolation of the velocity profile at the rotor and the 
rotor velocity imposed by the rheometer. Note that the 
bulk sample is sheared in all the cases reported here, 
proving that we only consider stresses larger than the 
yield stress in contrast with previous work [HES]. Lin¬ 
ear fits of the velocity profiles over 100 /im at a dis¬ 
tance of 50 /rm from the walls allow us to estimate the 
fluid velocity at the rotor and the stator, respectively 

fr) fs") 

noted Vj: ^ and uj, ^. Slip velocities are then computed as 
I, with i = r or s, and where Vw'^ = flR 

(s) 

and Vw = 0 , with fl the rotational speed of the rotor. 

(r) fs) 

The dependence of Vs and Vs ’ on the stress at the 
wall is reported in Fig. |^a). Slip velocities are plotted 
as a function of the viscous stress, a — ac, where cTc is 
fixed to the value of the yield stress inferred from the 
Herschel-Bulkley fit of the steady-state flow curve shown 
in Fig. I^c). Moreover, the local stress in the Taylor- 
Couette geometry is given by a(r) = T/[2TrH{R + r)^], 
where F is the torque on the rotor, so that in our ge¬ 
ometry the stress decreases by about 2e/R ~ 8 % from 
the rotor to th^tator. Such a variation is taken into 
account in Fig. by plotting and as a function 
of = cr(0) and = a{e) respectively. As a first key 
observation both slip velocities increase as a power law 
of the viscous stress over two decades of stress and four 
decades of velocities, with an exponent p = 2.17 ± 0.08. 
This result strikingly contrasts with previous analyses in 
which Vs was reported as a function of a rather than 
a — ac [m [El [23M5] . We have repeated these measure¬ 
ments at different temperatures above, around and below 
the jamming transition, respectively at T = 20 , 26 and 
34°C [Fig|^b-d)]. For all temperatures Vs robustly scales 
as a power law of the viscous stress (taking CTc = 0 for the 
liquidlike suspension at T = 34°C). Furthermore, the ex¬ 
ponent p is found to be the same at the rotor and at the 
stator within error bars, and to decrease smoothly from 2 
to 1 for increasing temperatures, taking the intermediate 
value p = 1.54±0.08 for T = 26°C [Fig.j^e)]. The deter¬ 
mination of p being sensitive to the value of a^ we have 
performed the following test to confirm the relevance of 
the above power-law scaling: instead of using the yield 
stress extracted from the Herschel-Bulkley fit of the flow 
curve, CTc is chosen as the value that minimizes the of 
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FIG. 4: (color online) Steady state slip velocity Vs vs the vis¬ 
cous stress a — ac for (a) a concentrated emulsion (0 ~ 0.75) 
sheared in a 3 mm gap Taylor Couette cell (data from [11] ') 
and (b) a 1% wt. carbopol microgel sheared in a rough mi- 
crochannel (data from [23]). In the absence of independent 
reliable estimates for ac, we used the yield stress obtained by 
minimizing the of linear fits of Inu^ vs ln((T — CTc), respec¬ 
tively ac = 19.7 Pa in the emulsion and ac = 5.5 Pa in the 
microgel. The red lines are the best power-law fits leading to 
p — 1.81 ± 0.06 and p — 1.63 ± 0.05 respectively. Insets in (a) 
and (b) show the exponents p vs the yield stress CTc inferred 
from a hundred least-square procedures on the same data sets 
to which a uniform noise of 10 % amplitude has been added. 
The dashed lines indicate the pair {ac,p) obtained in the ab¬ 
sence of noise. Colors code for the probability distribution of 
the yield stress values. 

linear fits of InUg vs ln((T — ac) when Uc is varied. The 
corresponding estimates are in good agreement with the 
previous values of p [Fig.j^e)]. Furthermore the values of 
CTc inferred from the least-square method are compatible 
within error bars with those from the Herschel-Bulkley 
fit (Table |T]). 

IV. DISCUSSION AND CONCLUSION 

We have shown that the power-law scaling Vg oc (a — 
ac)^ robustly describes the slip behavior of jammed sus¬ 
pensions of thermoresponsive pNIPAm particles sheared 
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above CTc- This prompts us to revisit two data sets where 
Vs was originally reported to scale as cr^: (i) a concen¬ 
trated emulsion sheared in a Taylor-Couette cell with 
smooth boundary conditions la and (ii) a carbopol mi¬ 
crogel sheared in a microchannel with rough boundary 
conditions [23]. Figure shows that both data sets are 
very well described by power laws of the viscous stress 
with p = 1.81 ± 0.06 (emulsion) and p = 1.63 ± 0.05 
(microgel). To further discriminate between Vg being a 
function of the stress or of the viscous stress, we de¬ 
termine the pairs of best fit parameters {<Jc,p) from a 
hundred realizations of the minimization procedure 
where a uniform noise of 10 % amplitude has been added 
to the data (insets of Fig. |^. In both cases the pair 
(tJc = 0 ,p ~ 2 ) lies at the very edge of the distribu¬ 
tion, which allows us to rule out a quadratic scaling with 
a. Our work thus provides an explanation to reconcile 
previous inconsistent results for stresses larger than the 
yield stress. It also draws a link with the slip behaviour 
of soft jammed materials below Uc where the material is 
solidlike and flows as a plug. Indeed, in microgels and 
emulsions for a < cTc, Vg was shown to increase either 
linearly or quadratically, i.e. Ug oc (tr — OgY■, where Og is 
the stress below which the material sticks to the bound¬ 
aries, with p = 1 or 2 depending on the nature of the 
interactions with the wall [nnn [331 [ssj. Our results 
show that for a > Og the apparent yield stress Og is re¬ 
placed by the bulk yield stress CTc and that the exponent 
p depends on the packing fraction. The latter result can 
be understood as follows: at low temperatures compared 
to the LOST, the exponent p = 2 can be interpreted as 
the result of a force balance on compressed particles ex¬ 
periencing Hertz-like deformation and lubrication forces 
with the walls, as first proposed for plug flows in mm 
At intermediate temperatures closer to the LOST, such 
a framework is probably no longer valid as the microgels 
may experience not only lubrication forces but also coop¬ 
erative rearrangements thanks to both the lower volume 
fraction and the non-zero local shear rate. Let us indeed 
recall that the material is sheared in the bulk for all data 
reported here, whatever the shear rate or the tempera¬ 
ture. The lesser importance of the lubrication interac¬ 
tions for decreasing packing fractions could thus account 
for the decrease of the exponent p, and for the transition 
to dilute-like behavior of the suspension, although the 
packing fraction is still above jamming. 

To conclude, we have shown that slip velocities in 


yield stress fluids made of soft jammed particles and 
sheared in the bulk increase as power laws of the viscous 
stress over a range of temperatures that encompasses 
both sides of the jamming transition. Therefore Vg can 
be written under the single form Ug oc [cr — a 
with ac — 0 and p ~ 1 below jamming and where 
both CTc and p are increasing functions of (j) above 
jamming. These results hint at a link between the bulk 
behavior of the material and the slip at the wall, as 
also evidenced by some of the nonlinear velocity profiles 
reported in Fig. Whether such nonlinearity should 
be attributed to some flow-concentration coupling m 
or to long-range effects of the boundaries associated 
with non-local effects [l2l [38l [39] remains an open issue. 
In any case a complete description of the present data 
shall certainly involve some coupling between bulk and 
boundary effects. Such a point of view contrasts with 
the interpretation of wall slip in foams, where slippage is 
mainly controlled by the physics at the boundaries and 
decoupled from bulk properties [TUI HOI HT] . Clearly a 
more detailed study of boundary-driven vs bulk-driven 
effects involving wall slip remains to be conducted by 
collecting data on various soft glassy materials including 
foams. Finally, our work paves the way for a more 
realistic modeling of flows of dense suspensions in the 
presence of slip. Spatially resolved model such as SGR 
and fluidity models [13] should focus on the scaling of Vg 
with a before turning to transient regimes for which wall 
slip goes hand in hand with long-lived heterogeneous 
dynamics [331H3] . 
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